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FOREWORD 

This  report  describes  the  results  of  an  inhouse  investigation 
dealing  with  an  analysis  of  the  response  of  aircraft  encountering 
aircraft  wake  turbulence.  The  work  was  partially  sponsored  by  The 
Pennsylvania  State  University  and  later  by  the  Air  Force  Flight 
Dynamics  Laboratory.  Beth  organizations  supplied  the  author  with 
computer  time. 

Special  recognition  should  be  given  to  Dis.  B.  W.  McCormick 
and  J.  J.  Eisenhuth  cf  The  Pennsylvania  State  University  and 
Dr.  G.  Kurylowich  of  the  Air  Force  Flight  Dynamics  Laboratory  for 
their  technical  guidance. 

This  investigation  was  performed  during  the  period  from 
September  1972  to  April  of  1974.  The  manuscript  was  submitted  by 
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SECTION  I 
INTRODUCTION 

The  wing  of  an  airplane  generates  from  its  trailing  edge  a  vortex 
sheet  which  rapidly  rolls  up  into  two  trailing  tip  vortices.  These 
vortices  nave  been  observed  to  persist  for  distances  of  up  to  10  miles 
or  more  behind  large  aircraft.  It  can  be  easily  shown  that  the 
strength  of  the  trailing  vortices  increases  with  the  size  of  the  air¬ 
craft  that  generates  them,  whereas  the  maximum  tangential  velocity 
increases  inversely  witn  the  velocity  of  the  aircraft.  Hence  the  most 
severe  vortex  system  to  be  expected  is  that  generated  by  a  large  air¬ 
craft  during  either  take-off  or  landing.  Therefore  the  vortex  hazard 
is  greatest  in  the  terminal  area.  The  probability  of  a  vortex 
related  accident  is  becoming  more  severe  with  the  growing  disparity  in 
aircraft  size  and  with  increasing  air  traffic.  During  the  past 
several  years,  there  has  been  a  renewed  interest  in  the  hazard  due  to 
aircraft  trailing  vortices.  Much  of  the  analytical  and  experimental 
research  conducted  during  this  period  deals  with  methods  of  predicting 
the  vortex  velocity  field  as  well  as  its  rate  of  decay.  There  are, 
however,  two  areas  which  have  not  received  much  attention.  These  are 
the  effect  of  atmospheric  conditions  on  the  transport  and  decay  cf  the 
vortices  and  the  dynamic  behavior  of  an  aircraft  penetrating  the 
vortex  system.  It  is  the  latter  topic  which  is  the  subject  of  this 
report. 

The  need  to  knew  the  dynamic  behavior  of  an  aircraft  encountering 
a  vortex  wake  becomes  clear  when  one  reviews  the  present  goals  of  the 
Federal  Aviation  Administration  (FAA).  The  FAA  plans  to  increase 
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airport  and  airway  capacity  by  a  factor  of  two  by  1980  and  by  a  factor 
of  five  by  1995.  In  order  to  accomplish  a  two  foxd  increase  by  iQ30 
the  FAA  oroposes  "o  improve  landing  aids,  reduce  separation  disc mces 
between  aircraft  to  two  miles,  and  to  utilize  dual  runways  vhicn  ar^ 
less  than  2500  feet  apart. ^  This  of  course  is  in  direct  oppos .  on 
to  the  suggestion  that  aircraft  be  separated  by  greater  aiscar  .  to 
insure  adequate  safety  from  aircraft  wake  turbulence.  The  FAA  .-cog¬ 
nizes  that  the  vortex  hazard  poses  a  serious  threat  to  aircraft 
safety  and  airport  utilization  and  is  attempting  to  find  a  sole'  -„en 
b>  supporting  research  in  toe  areas  of  vortex  dissipation  and  vertex 
detection. 

Efforts  to  eliminate  the  vertex  hazard  by  causing  the  vortices 
tc  dfceay  earlier  have  not  as  yet  been  successful.  The  various  con¬ 
cepts  which  have  been  considered  include  wing  ues-igr.,  mass  injection, 
spoilers,  oscillating  control  surfaces,  engine  location,  and  trailing 
bodies.  The  subject  of  vortex  dissipatocs  is  treated  in  reference  2. 
All  of  these  devices  change  the  velocity  field  in  the  vicinity  of  the 
wing.  However,  these  devices  do  not  appear  to  change  the  vortex 
system  at  distances  far  downstream  of  the  aircraft.  For  example,  it 
was  suggested  by  Corsiglia,  Jacobsen,  and  Cbigier*  that  placing  a  small 
vertical  panel  near  the  wing  tip  would  cause  a  modification  of  the 
vortex.  Wind  tunnel  results  seemed  to  verify  their  idea.  However, 
when  the  device  was  placed  on  a  Convair  990  the  results  were  incon¬ 
clusive.  The  pilots  of  the  aircraft  probing  the  990’s  wake  reported 
that  they  were  unable  to  detect  any  difference  between  the  Cor.v.tir  990 
with  or  without  the  vortex  dissipator. 
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The.  leading  candidate  for  a  vortex  detect* ta  system  is  an  acoustic 
sensor.  There  are  several  devices  presently  under  consideration.  They 
all  have  limited  range  and  can  only  determine  the  vortices  located  in 
a  very  narrow  field  of  view.  Figure  1  taken  from  reference  1  shows  an 
artist’s  conception  of  the  planned  vortex  detection  system. 

Even  if  an  operational  vortex  detection  system  is  developed 
within  the  next  few  years  one  still  needs  to  know  how  hazardous  a 
detected  vortex  is  to  the  next  aircraft  to  pass  through  the  same  air¬ 
space.  The  information  included  in  this  dissertation  will  aid  in 
determining  safe  operating  procedures  in  »"he  terminal  area. 

This  report  attempts  to  evaluate  the  reiai ionship  between 
the  motion  of  the  penetrating  aircraft  and  the  size  and  operating  con¬ 
ditions  of  the  generating  aircraft.  It  also  considers  the  age  of  the 
vortex  and  the  manner  in  which  the  vortex  is  encountered. 
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SECTION  II 

PREVIOUS  INVESTIGATIONS 

Tv.'1,  problem  associated  with  an  aircraft  penetrating  t'r.e  trailing 
vortex  system  of  another  aircraft  was  first  analyzed  by  Bleviss  In 
his  analysis  he  determined  the  aaxiss'm  roll  rate  to  be  expected  if  the 
aircraft  were  suddenly  immersed  along  the  axis  of  one  of  the  vortices 
as  shown  in  rigure  2.  The  induced  roll  rate  was  compare.-’  with  the 
roll  rate  available  using  aileron  control.  The  analysis  pointed  out 
the  hazard  to  light  aircraft  encount  .ring  the  wake  of  a  civilian 
transport  aircraft. 

Introduction  of  let  transports  into  commercial  aviation  in  the 

late  fifties  prompted  a  renewed  interest  in  the  vortex  hazard. 

4,  5,  6  7 

Several  papers  by  McGowan  and  Vetiaore  and  Reeder  were  pub- 

published  in  the  sixties.  McGowan  analyzed  the  behavior  of  ar.  air¬ 
craft  traversing  the  vortex  system  (see  Figure  2).  He  was  the  first 
to  include  unsteady  aerodynamics  into  the  modeling  of  the  problem. 

Some  typical  results  from  his  analysis  are  shown  in  Figure  3.  His 
calculations  indicated  that  light  normal-  category  and  light  trans¬ 
port  -  category  aircraft  could  experience  loading  conditions  that 
exceed  the  design  limit  and  in  some  cases  the  design  ultimate  load 
factors.  He  also  concluded  that  the  pilot  would  aggravate  the 
loading  condition  by  his  phasing  of  the  elevator  deflection. 

Wetocre  and  Reeder  used  HcGownn’s  and  Bleviss*  results  to 
develop  procedures  in  the.  airport  terminal  area.  Tneir  main  con¬ 
clusion  was  that  the  vortex  hazard  could  be  substantially  reduced  by 
developing  air-traffic  conrrol  procedures  which  emphasized  appropriate 
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spacing  ami  control  of  flight  paths. 

In  the  late  sixties  and  early  seventies  the  interest  in  vortex 
wakes  was  renewed  with  the  advent  of  the  so  called  "Jumbo  Jet". 
Several  symposiums  were  held  to  express  ideas  on  how  to  solve  the 
vortex  problem.  Most  of  the  napers  presented  at  these  svnposia 
dealt  with  the  theories  of  vortex  decay,  atmospheric  effect  or  the 
vortices  and  means  of  alleviating  the  vertex  hazard  by  hastening  he 
vortex  breakup.  Several  papers  on  the  dynamics  of  a  vortex  encounter 

were  presented  by  Hackett  and  Theisen^,  Andrews*' ^  and  Condit  and 

_  11 
Tracy  . 

Hackett  and  Theisen  presented  results  from  a  computer  simulation 
of  a  vortex  encounter.  Their  results  indicated  the  significant 
effect  that  the  vertical  tail  has  on  the  vortex  upset.  Figure  4 
taken  from  their  paper  shows  the  response  of  a  Jet  transport  tc  • a 
encounter  with  the  wake  of  a  Jumbo  Jet.  Their  results  also  indicate 
that  the  pilot's  reaction  could  aggravate  the  vortex  upset.  This  can 
be  explained  by  examining  Figure  4.  As  the  alrcraft  enters  the 
vortex  from  the  top  side  the  aircraft  starts  to  roll  in  a  counter 
clockwise  direction.  The  vertical  tail  is  subjected  to  a  lateral 
velocity  ^  the  left  causing  the  aircraft  to  yaw  to  the  right.  As 
the  aircraft  moves  to  the  right  the  left  wing  experiences  the  high 
velocity  upwash  of  the  right  vortex,  "his  causes  the  aircraft  to 
roll  in  a  clockwise  direction  which  is  opposite  of  the  initial  roll 
disturbance.  Thus  the  pilot's  attempt  to  control  the  aircraft  would 
be  momentarily  out  of  phase  with  the  motion  of  the  aircraft. 

Andrews  presented  the  results  of  tests  conducted  at  the  NASi 
Flight  Research  Center  in  which  probe  aircraft  were  flown  into  the 
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wakes  of  heavy  Jumbo  Jet-type  aircraft.  Test  results  indicated  that 
aircraft  with  short  wing  spans  could  experience  uncontrolled  upsets 
for  up  to  8  nautical  miles  behind  a  Jumbo  Jet.  They  also  concJ  Ued 
that  the  maximum  induced  rolling  moments  exce<ded  the  lateral  cmtrcl 
pcwfr  most  of  the  probe  aircraft*  when  they  were  wither.  the 
minimum  separation  distances  normally  maintained  during  take-.  and 
landing  operations. 

Condit  and  Tracy  presented  the  results  of  3oeing's  wake  ti-bu- 
lence  Investigations.  Based  on  their  tests  and  analyses  they 
recommended  that  the  IAA  use  a  5  mile  separation  for  iiiht  airplanes 
following  heavy  transports.  They  also  concluded  that  the  74?  and  707 
produced  similar  dynamic  responses  to  aircraft  encountering  their 
respective  wakes. 

The  most  recent  works  on  aircraft  vortex  interaction  are  by 
12  13 

Bernstein  and  Iversen  and  by  Nelson  and  McCormick  .  Beiustezr.  and 
Ivarsen  developed  a  3  degree  of  freedom  analog  simulation  of  the 
vortex  problem.  They  were  primarily  concerned  with  penetrations 
along  the  axis  of  the  vortex  cores.  Their  results  showed  that  the 
response  of  a  C-130  encountering  the  wake  of  a  C-5A  at  a  distance  of 
2  nautical  miles  was  more  severe  than  the  response  of  a  C-130 
260  feet  behind  another  C-130.  Thus  their  results  indicate  the 
importance  of  the  relative  size  of  the  penetrating  aircraft  to  that 

*  The  probe  aircraft  were  a  Convair  990,  DC-9,  Lear jet  and  a 
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of  the  generating  aircraft. 

Nelson  and  McCormick  shewed  that  relatively  large  aircraft  can 
be  susceptible  to  vortices  generated  by  large  jet  transports.  Their 
conclusions  were  based  upon  a  review  of  accident  records  as  well  as  a 
computer  simulation  of  the  aircraft  vortex  penetration.  The  com¬ 
puter  simulation  consisted  of  the  equations  of  motion  with  6  degrees 
of  freedom  and  also  included  control  input  Dy  the  pilot.  Figure  5 
taken  from  reference  13  shows  the  response  of  a  light  jet  transport 
(DC-9)  to  an  encounter  with  the  wake  of  a  Jumbo  Jet  (DC-10).  It 
should  be  noted  that  the  maximum  roll  argle  exceeded  50  degrees  even 
with  pilot  feedback. 

This  brief  survey  of  the  papers  dealing  with  aircraft  vortex 
interaction  is  by  no  means  complete.  There  are  other  papers  pub¬ 
lished  on  this  topic  and  the  reader  is  referred  to  a  paper  by 
14 

McCormick  for  an  extensive  bibliography  on  trailing  vortex  sys¬ 
tems  and  aircraft  responses.  The  papers  presented  in  this  section  are 
in  the  author's  opinion  the  most  significant  ones  dealing  with  the 
dynamics  of  a  vortex  penetration. 
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SECTION  III 
ANALYTICAL  STUDY 

To  investigate  the  dynamic  behavior  of  an  aircraft  penetrating 
a  vortex  wake  a  computer  simulation  wa&  developed.  The  program  was 
written  in  the  FORTRAN  IV  language  and  consists  of  nine  subroutines 
and  seven  function  subprograms.  In  order  to  accurately  model  the 
penetration  problem  the  computer  program  had  to  include  the 
following: 

i  a  realistic  description  of  the  vortex  velocity  field 

ii  a  means  of  determining  the  induced  velocity  distribution 
affecting  the  penetrating  aircraft 

iii  a  means  of  computing  the  unsteady  aerodynamics  associated 
with  the  vortex  velocity  field 

iv  the  complete  set  of  rigid  body  equations  of  motion 

v  and  finally  a  means  of  introducing  pilot  control. 

Each  of  the  above  topics  is  discussed  in  detail  In  the  following 
sections.  Included  in  the  discussion  are  the  reasons  for  the  selec¬ 
tion  of  one  theory  over  another  as  well  as  the  assumptions  made  in 
the  equations  used  in  the  program. 

Vortex  Decay 

The  rate  of  decay  of  the  velocities  within  the  vortex  wake  at 
a  given  distance  behind  the  generating  aircraft  are  difficult  to 
calculate  and  to  a  certain  extent  largely  unknown.  There  are 
numerous  theories  which  c/n  be  used  to  calculate  the  velocity 
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distribution.  However,  lacst  of  these  theories  prove  to  be  inadequate 
when  compared  with  experimental  data.  Several  theories  for  vortex 
decay  will  be  discussed  and  some  of  the  theories  will  be  compared 
with  experimental  data  obtained  by  the  Federal  Aviation  Administra¬ 
tion  at  their  National  Aviation  Facilities  Experimental  Center  (NAFEC) . 

The  most  familiar  solution  for  the  diffusion  of  a  line  vortex  is 
that  obtained  by  Lamb^. 


V„  (r , t) 

o 


] 


where  ?  is  the  circulation  and  v  is  the  kinematic  viscosity.  Various 

experimenters  such  as  Rose  and  Dee^^,  Bisgood,  Maltby  and  Dee^,  and 
18 

Squire  have  tried  to  use  a  modified  form  of  Lamb's  solution, 

.  _ r! 

r  4(\*c)t 

V?<r’t)  '  2Tr  [1-c  1 

where  c  is  the  eddy  viscosity.  The  value  of  the  eddy  viscosity  is 
assumed  to  be  proportional  to  the  circulation,  i.e.. 


z  •  a: 

The  constant  of  proportionality,  however,  has  proved  to  be  very 

-a 

difficult  to  measure  and  in  fact  has  been  found  to  lie  between  10 
and  10  .  On  the  other  hand  Newman  and  Dosanjh  u  replaced  the  kine¬ 

matic  viscosity  by  an  eddy  viscosity  v^.  Newman  found  that  the  ratio 
of  eddy  viscosity  to  the  kinematic  viscosity  had  to  be  varied  from 

one  axial  station  to  another  in  order  to  match  the  data.  Dosa.ijh  found 

good  correlai _on  with  his  data  by  using  an  eddy  viscosity  of  *=  lOv. 
21 

Recently,  Rurylovich  applied  the  •’ame  type  of  analysis  to  flight 

test  data  and  arrived  at  a  value  for  vx  of  260v/cos^\  ,,  where  X 

i  c/4  c/4 
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chord  line.  Kurylowich's 
radius,  and  tangential  velo- 


r  <  2a 


JThf  additional  terms  ir.  the  preceding  expressions  are  included  to  take 

into  account  the  effect  of  mass  flow  injection  into  the  vortex  core. 

Inclusion  of  these  terms  was  based  upon  results  obtained  by 
22 

Poppleton  .  Poppleton  found  that  the  core  radius  and  turbulence 

levels  increased  as  the  mass  flow  injected  into  the  vortex  core  va 

increased.  Therefore  the  constants ,  k£  and  k^,  which  must  be 

uetermined  experimentally  are  included  to  take  into  account  the 

increased  eddy  viscosity  due  to  mass  injection. 

23 

McCormick,  Tangier,  and  Sherrieb  showed  chat  the  trailing 
vortices  could  be  described  by  geometric  similarity  considerations. 
Their  work  was  based  upon  analysis  of  wind  tunnel  and  flight  test 
data.  The  express ons  developed  in  their  analysis  are 

Va  (a)  -  .53  CL  ” 

'  o 

F(a)  =  •  Ifcl'a- 


is  the  sweep  angle  of  the  wing  quarter 

equations  for  the  eddy  viscosity,  core 

city  are  as  follows: 

260v  ,,  ,,  , .  2^ 
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1  +  ln(r/a) 
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1  +  In  r/a 
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where 

ve0(a) 
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Ve(r) 

r(a) 

r 

00 
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the  maximum  tangential  velocity  immediately  behind 
the  wing. 

the  maximum  tangential  velocity  at  a  distance  z 

downstream  of  the  aircraft. 

the  tangential  velocity  at  any  radius, 

the  circulation  at  the  core  radius. 

the  total  circulation. 

the  circulation  3t  any  radius. 


The  expressions  for  Vg(a)/Vgo(a)  and  T(a)  were  determined  from  wind 

tunnel  and  flight  test  data.  However,  the  expression  used  for  the 

24 

circulation  is  based  upon  the  analysis  of  Hoffman  and  Joubert 
Their  analysis  predicts  a  logarithmic  variation  of  the  circulation 
with  the  radius. 

In  a  later  paper  dealing  with  an  analysis  of  experimental 

measurements  of  the  vortex  wake  behind  large  jet  transports  Eisenhuth, 

25 

McCormick,  Nelson  and  Garode  suggested  that  T(a )/F  depends  upon  the 
lift  coefficient  instead  of  being  a  constant.  Figure  6  shows  a  plot 
of  TCaj/r^  versus  which  does  indeed  appear  to  varv  in  a  linear 
manner.  Eise.ihuth,  et  al,  suggest  that  the  following  equations  should 
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be  used  to  predict  the  velocity  distribution  through  the  vortex. 


r(a) 


.17 


1  /? 

—5 —  -  .02(1  +  .00063  z/c  CL)  ' 
coCi0 

J  C£ 

r (a)  -  .085  Cj/  c0(-~)  V 


V0(a) 


.63  CL  V  (1  +  .00063 


-1/2 


Ve(r)  ^  1  4-  In  r/i 
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In  several  recent  papers  Donaldson^  *  has  proposed  using 
28 

Betz's  analysis  to  describe  the  roll  up  process  as  well  as  the 
velocity  distribution  through  the  vortex  system.  Betz  used  the 
following  facts  in  his  analysis. 

1.  All  the  vorticity  shed  by  each  half  of  the  wing  is 
found  rolled  up  into  the  trailing  vortex  behind 
the  appropriate  half  of  the  wing. 

2.  The  "center  of  gravity"  of  the  vorticity  distribu¬ 
tion  remains  at  a  constant  distance  from  tne  plane 
of  lateral  symmetry. 

3.  The  "moment  of  inertia"  of  the  vorticity  shed  by  each 
half  of  the  wing  about  its  "center  of  gravity"  is 

a  constant. 

Rewriting  these  statements  in  mathematical  form  fields. 
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(for  elliptic  distribution) 
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Betz  used  these  expressions  to  determine  a  circulation  distribution  as 
shown  in  Figure  7.  Donaldson  developed  an  approximation  to  this  dis¬ 
tribution  which  is 


~  -  [6 (r/b)  -  9(r/b)V/2 

1  o 


0  <  r/b  <  1/3 


~  =  1  r/b  >  1/3 
1  o 


Stated  in  terms  of  the  tangential  velocity  the  expressions  are 


Vfi  «  [6(r/b)  -  9 (r/b ) 2 ] 1/2  0  <  r/b  <  1/3 

r0 

V  -  r/b  >  1/3 

zirr 

In  a  later  paper  Donaldson  extended  the  preceding  analysis  to  include 

aircraft  having  flaps  deflected.  Figure  8  taken  from  reference  27 

shows  a  comparison  of  his  technique  with  experimental  data. 

The  methods  discussed  in  this  section  were  compared  with  experi- 

20  in  12 

■  ital  data  obtained  by  the  FAA's  Wake  Turbulence  program.  *  *  ’ 

Data  was  obtained  by  flying  aircraft  of  various  sizes  past  an  instru- 
-  "ted  tower.  The  tower  was  equipped  with  hot  wire  and  hot  film 
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velocity  sensors.  Vortices  of  different  ages  were  obtained  by 
appropriately  spacing  the  aircraft  vertically  and  laterally  so  that 
it  would  tare  a  longer  time  for  the  vortices  to  reach  the  tower. 

Figure  10  shows  a  sketch  of  the  tower  fly-by  technique.  It  should  be 
mentioned  that  the  data  was  accumulated  in  ground  effect.  This  was 
due  to  the  relatively  low  height  of  the  test  tower.  Fui  more  informa¬ 
tion  on  the  test  facilities  the  reader  is  referred  to  a  paper  pre- 
33 

sented  by  Garodz  which  describes  the  test  procedures  in  detail. 

Figures  10  a-h  show  the  comparison  between  the  vortex  theories 

and  the  experimental  data.  The  airplanes  used  to  generate  these  data 

range  from  a  DC-6B  up  to  a  wide  body  jet  tiansport.  The  vortex  ages 

range  from  7.7  to  31  seconds  which  correspond  to  distances  of  1770  to 

7150  feet  behind  the  generating  aircraft.  Only  the  theories  which 

compared  favorably  with  the  data  were  plotted  on  the  figures. 

25 

Examination  of  the  figures  show  that  Eisenhuth's  expressions  consis¬ 
tently  match  the  data  with  the  exception  of  the  747  runs.  It  should 

25 

be  recalled  that  the  expressions  proposed  by  Eisenhuth,  et  al  were 
based  on  McCormick's"^  analysis.  The  disagreement  between  the  pre¬ 
dicted  and  measured  data  for  the  747  is  not  easily  explained.  It  has 
been  suggested  in  the  literature  that  the  exhaust  from  the  ?47's  out¬ 
board  engines  is  entrained  in  the  vortex  and  alters  its  structure. 

This  seems  plausible  when  one  examines  the  data  obtained  by  Marchman 
34  22 

and  Mason  or  Poppleton  .  These  investigators  studied  the  effect 
of  air  injection  into  the  vortex  core.  Figures  11  and  12  taken  from 
reference  34  show  the  effect  of  axial  blowing  (injection)  on  the  velo¬ 
city  distribution,  vortex  core,  and  circulation.  The  data  show  that  blowing 
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into  the  vortex  core  reduces  the  maximum  tangential  velocity  and 
increases  the  core  radius.  The  circulation  distribution  is  altered 
in  the  vicinity  of  the  vortex  core,  however*,  at  radial  distances 
greater  than  the  vortex  radius,  the  circulations  approach  the  theoreti¬ 
cal  bound  vortex  strength.  Also  included  in  Figure  11  is  a  prediction 
based  on  Kurylowich's  method.  It  should  be  no*-ed  that  the  eddy 
viscosity  had  to  be  adjusted  to  a  value  of  20v.  These  figures  were 
included  to  help  explain  the  difference  between  the  data  of  the  777 
and  the  other  aircraft.  In  a  later  section  the  effect  of  modifying 
the  vortex  velocity  field  on  the  dynamic  response  of  a  following  air¬ 
craft  will  be  discussed. 

Based  on  the  comparisons  presented  and  others  not  included  in  this 
report  it  was  concluded  that  Eisenhuth's  expressions  could  be  used  to 
accurately  model  the  velocity  distribution  through  the  vortex  system. 

Vortex  Induced  Velocities 

The  wing  is  divided  into  N  strips  as  shown  in  Figure  13.  The 
velocity  induced  by  the  trailing  vortex  system  is  resolved  into  com¬ 
ponents  [U  (t),  V  (t),  W  ( t) ]  at  control  puints  located  at  the  j/4 
g  g  g 

chord  point  of  each  strip  on  the  wing  and  tail  surfaces. 

The  procedure  for  resolving  the  induced  velocity  into  its  compo¬ 
nent  .  [U  (t),  V  (t),  V  (c) j  is  as  follows: 
g  g  S 

1.  The  inertial  axes  system  is  aligned  so  that  the  X  axis  is 
parallel  to  the  axis  of  the  vortex  filaments.  (See  Figure  14) 

2.  Using  the  transformation  shown  below  the  control  points  c3n 


be  resolved  into  the  inertial  frame. 


Vortex  system  parallel  to 
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3.  The  perpendicular  distance  of  each  control  point  fro-n  the 
vortex  filaments  is  equal  to  the  following  expressions. 

YJ  =  Ycg  +  Ycpi 


Zj  -  Zcg  +  ZCpj 

Rii  -  ((Yi  -  Yv)2  +  (ZA  -  Zv)2} 


2,1/2 


R2j  “  UYi  -  YVd)2  +  (Zi  -  Zv)2}1/2 

4.  The  velocity  components  in  the  iuertial  frame  are 
VYv  Yi-Yvd 

■  X  (V  - ^  (^r> 


■  -VelJ  ^  +  V*2j  ^ 

if  the  control  point  lies  outside  the  core,  and 
Y4-Yv 

“8,  *  V*lj 

v  -  -ve  (fi±) 

»j  lj  a 


port  vortex  core 


—  .Y1-Yvd 

V  *  W  W~> 


z1_zv 

V  "  V02j  (J7-} 


starboard  vortex  core 
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if  the  control  point  lies  within  the  vortex  core. 

5.  The  velocity  components  referred  to  the  body  axes  system  are 
obtained  by  using  the  following  transformation. 


_ 

" 

ugj 

re  c  c  -c 

~0  $  ‘6 

V^, 

at 

S,S„C  S  8  C, 
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V  $  V  $ 

J 

However,  Ug^  is  assumed  to  be  zero,  i.e.,  zero  axial  '  elocity  in 
vortex  field.  This  allows  the  expressions  to  be  reduced  to  the 
following  form. 
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Aerodynamics 

The  problem  of  determining  the  aerodynamic  forces  acting  on  an 
aircraft  in  a  time-varying  velocity  field  is  indeed  quite  formidable. 
However,  in  order  to  accurately  predict  the  response  of  an  aircraft 
penetrating  a  vortex  system,  some  means  nf  incorporating  the  unsteady 
aerodynamics  into  the  equations  of  motion  of  the  aircraft  is  neces¬ 


sary. 
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There  have  of  course  beer,  numerous  lifting  surface  theories 
developed  over  the  years.  I-  the  early  sixties  two  techniques 
dominate  the  litr  tur._ ,  they  are  the  voitex  and  doublet  lattice 
methods.  u.e  use  of  discrete  vortices  for  the  solution  of  the 
s,eady  lifting  surface  problem  was  first  used  by  Falkner^.  Falkner, 
however,  had  to  make  assumptions  on  the  load  attributions  in  order 
to  keep  the  number  of  unknowns  to  a  minimum.  In  the  sixties, 


"*6  37  38  3°  uu 

Hednan  ,  Giesin*  ,  Rubbert  .Woodward  Alb  an  c  and  Rodden  and 

others  developed  lifting  surface  theories  which  did  not  make  any 

prior  assumptions  on  the  load  distribution.  This  was  possible  due 

to  the  feasibility'  of  using  modern  digital  computers  to  invert  the 

matrices  associated  with  the  finite  element  methods. 

The  method  discussed  in  this  section  is  an  unsteady  vortex 

41 

lattice  technique.  This  tecnnique  is  an  extension  of  Rudhman's 
work  to  the  finite  wing. 

Now  consider  the  integral  equation  of  the  linearized  lifting 

42 

surface  theory  as  formulated  by  Reissner  and  others, 
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t)  [x-;]+iw(4,-,  t)  [y—  ] 
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wnere  Rd  -  wiug  region 
*  wake  region 

Y  &  c  -  the  circulation  per  unit  chord  and  per  unit  an, 
r ■  :~'t ivelv. 
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See  Figure  15  for  clarification  of  the  above  symbols.  The  integral 
equation  relates  the  downwash  distribution  to  the  unknown  running 
circulations.  The  boundary  condition  which  mu3t  be  satisfied  is  as 


3Z 

c.y.o.t)  •»  (x,y,o,t)  +  U 

3t 


a^a^»y  i 
3x 


This  is  a  mathematical  statement  of  the  flow  tar»ency  condition. 

If  a  discrete  vortex  method  is  used  the  above  integral  equation 
is  replaced  by  a  system  of  simultaneous  algebraic  equations. 

wi  *  ^  r3 

or  in  matrix  notation 
{w}  -  [A]{r> 

where 

{W}  -  {Wg}  +  {*„} 

[”3}  *  normalwash  which  is  known  from  the  boundary  condition. 
{Ww}  “  upwash  (downwash)  due  to  vortices  shed  into  the  wake. 

[A]  °  matrix  of  aerodynamic  influence  coefficients  relating 

the  upwash  (downwash)  at  a  point  i  due  to  a  singularity 
at  point  j. 

{r}  «  ..  vseshoe  vortex  strength. 

In  the  above  formulation  both  {Wg}  and  {Ww }  are  known.  {Wg}  is 
determined  from  the  boundary  condition.  {Ww }  is  determined  from  the 
following  equation 

-  l  [B]fc  {r}k 

k=l 


where 


[B]k  =  matrix  of  Influence  coefficients  relating  the  upwash 
(downwash)  at  a  point  i  (on  the  v  r  g  •  :uie  to  a 


Y(x,y,t)  and  f(x,y,t)  are  the 
circulations  per  unit  length  in 
the  x  and  y  directions 


Figure  15.  Vortex  Model  of  the  Finite  Wing 
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singularity  at  point  j  (in  the  wake). 

{F)k  =  strength  of  vortices  shed  into  the  .ake. 
n  =  number  of  time  steps  or  {T}'s  shed  into  the  wake. 

It  should  be  noted  that  (?V  is  known.  This  is  so  because  the  circu¬ 
lation  shed  into  the  wake  is  related  to  the  cnange  in  -i-g  circula¬ 
tion.  That  is  the  wake  circulation  shed  during  a  time  interva'  t  is 
equal  and  opposite  to  the  change  in  wing  circulation  during  the  me 
time  interval. 

The  approach  used  in  solving  the  unsteady  lifting  surface  problem 
is  to  represent  the  wing  by  a  network  of  closed  horseshoe  vortices 
uho-ch  are  distributed  in  both  the  chordwise  and  spanwise  directions. 
The  wake  is  constructed  as  shown  by  the  simple  sketch  in  Figure  16. 

The  strength  of  the  individual  vortices  are  determined  by  applying 
the  wing  boundary  condition  at  as  many  control  points  on  the  wing  as 
there  are  vertices.  The  resulting  system  of  simultaneous  equations 
is  then  solved  to  determine  the  strength  of  the  vortex  filaments. 

The  spacing  and  location  of  vortex  filaments  is  shewn  in  Figure 

17.  Notice  that  the  bound  vortices  are  located  at  the  1/4  cnord  of 

each  chordwise  panel.  The  control  point  is  located  at  the  3/4  chord 

point  of  the  panels.  This  choice  of  control  point  location  has  been 
A3  A  A 

shown  by  James  and  De  Young  to  be  the  optimum  for  the  two  dimen¬ 
sional  cases.  Although  nothing  has  been  said  about  imposing  the 
Ku-rta  condition  in  this  method,  it  has  been  found  by  numerical 
experimentation  that  the  Kutta  condition  will  be  satisfied  when  the 
control  point  is  located  at  the  3/4  chord  location.  The  location  of 
the  trailing  filament  was  obtained  by  using  the  following  trans¬ 


formation. 
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Figure  16.  Vortex  Wake 


Vortex  -  1 / 4  Chord 

X 


Control  point  -  3/4  Chord 


Figure  17.  Chordvise  Vortex  Spacing 


Figure  18.  Vortex  Spacing  Spanvise 
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9V 


fVy  -  cos  6V 


Multhopp's  transformation 


where 

m  *  number  of  stations  (spanwise) 
v  *  identifying  subscript  of  a  particular  station 
nv  *  nondimens ional  distance  n  “  y/^/2 

For  further  clarification  see  Figure  18.  The  influence  coefficient 
matrices  are  determined  by  using  the  integrated  form  of  the  Biot- 
Savart  law  for  a  straight  line  segment.  Applying  this  equation  to 
each  side  of  the  closed  horseshoe  filament  one  can  determine  the 
induced  velocity  at  any  control  point  due  to  the  vortex  box. 

The  algorithm  used  to  solve  for  the  unknown  circulations  is 
listed  below. 

Step  1.  Filaments  are  placed  on  the  wing  as  previously 
described. 

Step  2.  Influence  coefficients  [A}  are  determined. 

Step  3.  {W}  is  computed.  For  example  the  normalvash  and  wake 

induced  normal  velocity. 

Step  4.  Unknown  circulations  are  determined 
{D  *  [A]_1{W] 

S Lop  5.  Shed  vortices  are  moved  downstream  by  the  distance  VAt. 
Step  6 .  Repeat  steps  3  through  5. 

The  above  numerical  procedure  yields  the  circulation  or.  the  wing  at 


discrete  points  in  time.  The  lift  on  the  wing  panels  is  determined 
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by  integrating  the  pressure  difference. 

a  r 


AP  “  pVy(x)  +  p 


at 


Y(x)  dx 


This  equation  is  obtained  from  the  unsteady  Bernoulli  equation. 
Expressing  Euler's  equation  in  victor  form 

||  +  I  ,,2  .  .  I  V  P 


and  the  velocity  in  teras  of  the  velocity  potential 
q  *  grad  6 

then  the  following  equation  can  be  obtained. 


V 

1  at 


+ 


ai 

2 


o 


Integrating  this  expression  yields 


at 


+ 


const. 


The  above  is  known  as  Kelvin's  equation  or  the  unsteady  Bernoulli 
equation.  Note  that  this  equation  only  applies  along  a  streamline. 
If  this  equation  is  applied  to  the  thin  airfoil  shown  in  Figure  21, 
one  can  determine  an  expression  for  the  unsteady  lift.  Consider  a 
small  segment  of  the  above  airfoil  section  shown  in  Figure  22. 

The  circulation  around  this  segment  is 
ydx  *  2AV  dx 
or 

AV  -  y(x)/2 

The  definition  of  the  potential  between  two  points  A  and  B  is  given 
by 
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2AVdx 

■SsL 

2 


Figure  22. 


Segment  of  Vortex  Sheet 
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4>(B)  -  ♦(A)  = 


Vs  ds 


Since  the  potential  at  the  leading  edge  of  the  airfoil  can  be 
arbitrarily  set  equal  to  zero  ($(L.E.)  =  0)  the  above  expression 
reduces  to 


r 

J  O 


4>  (x)  -  /  V8  ds 


or 


$«(*>  * 


*£(x) 


(V  +  xi2£>)  d* 


(V  -  3^  )  d. 


tfiere  u  and  £  indicate  the  upper  and  lower  surfaces  respectively. 
Substituting  back  Into  the  unsteady  Bernoulli  equation  yields  the 
pressure  difference. 


or 


pt  -  PU  -  -  »  jf  <♦»*♦»>  -  S/2  <V-vu2> 

r» 

A?  -  pVy(x)  -f-  p  ~  !  y(x)  dx 

«'  o 

The  lift  on  the  panel  can  now  be  found  by  integrating  the  pressure 
difference  from  the  leading  edge  to  the  trailing  edge. 


i  ^  g 

AP  dx  *  pV  /  y(x)  dx  +  p  — 


x 


£  -  pvr(t)  +  p 


?t 


v  J 
/' 


o  J  o 


y(x)  dx  dx 


T(x)  dx 
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The  equations  used  for  calculating  the  lift  for  a  particular  station 
are  shown  for  five  chordwise  vortices. 

„  ri  .  r1(t+  t)  -  ri(t) 

AP1  “  pV  Zxl  +  p  "At - 

r2  rx(t+  t)+r2(t+  t)-ri(t)-r2(t) 
ap2  -  ov  +  c  st 


ap< 


The  lift  per  unit  span  is  as  follows. 


Z  -  APjAxjl  +  + 


AP5 


x6 


£  -  APi(Ax1+^i)+  £P2  (^1  + 


Ax3  £x4  AX4  Axs 

+  AP3(  2  +  2  )  +  AP4O  2"  "+  2  ) 


AX5  Axa 

+  1P^"T 


■) 


The  moment  can  also  be  found  in  a  similar  manner.  The  total  lift  and 
moment  are  calculated  by  summing  the  increments  of  lift  in  the  span- 
wise  direction. 

The  technique  described  in  this  section  could  be  used  to 
determine  che  induced  drag.  However,  the  interpretation  of  the  drag 
results  computed  from  lifting  surface  theories  still  remains  a 
question  of  some  uncertainty  in  the  subsonic  speed  regime.  Another 
point  that  should  be  made  is  that  the  technique  does  not  take  into 


account  thickness  effects.  This  could  be  taken  into  account  by  usin'., 
source  distributions. 

The  unsteady  vortex  lattice  method  was  programmed  and  various 
test  cases  wehe  run.  Comparisons  were  made  with  both  two-  and  three- 
dimensional  results.  Figures  23  and  24  show  the  excellent  agreement 
between  the  classical  two-dimensional  solution  of  Theodorsen,  *u  -er 
and  Kussner  with  the  results  calculated  from  the  vortex  lattice  pro¬ 
gram.  There  is  a  difference  between  the  computed  values  and  the 
tneoretical  Kussner  and  Wagner  functions  for  the  first  few  semichord 
lengths  traveled  by  the  airfoil.  However,  the  Kussner  and  Wagner 
functions  describe  only  the  circulatory  lift  due  to  a  gust  penetra¬ 
tion  and  a  step  change  in  angle  of  attack  resoectively,  whereas,  the 
computations  based  on  the  equations  developed  for  the  vortex  lattice 
method  include  both  the  noncirculatorv  and  cir'ulatory  lift.  Since 
the  added  mass  effects  are  important  at  the  beginning  of  an  impulsive 
start  then  this  would  account  for  the  difference  between  the  curves. 
Figures  25  and  ?j  show  the  results  for  a  finite  wing  having  an  aspect 
ratio  of  6.  Again  the  difference  can  be  explained  by  t’.e  absence  of 
the  added  mass  effects  in  the  classical  solution. 

An  interesting  comparison  between  a  simple  strip  theory  and  the 
more  exact  lifting  surface  theory  is  presented  In  Figure  27.  This 
figure  shows  a  comparison  for  a  step  change  in  gust  velocity  (Kussner 
Problem).  The  difference  in  the  lift  coefficient  is  almost  negiigibi 
Figure  28  shows  a  comparison  of  the  roll  moment  coefficient  computed 
by  both  techniques  for  a  gust  field  which  /aries  linearly  from  the 
wing  centerline  to  rhe  wing  tip. 

Although  the  vortex  lattice  method  yields  extremely  accurate 


Theoret i 


.lBon  of  Unsteady  Vortex  Lattice  with  Kussner  and 
r  Functions 


Semirhords  traveled 
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results,  the  feasibility  of  using  a  complete  lifting  surface  theory 
to  calculate  the  loads  on  an  aircraft  penetrating  a  trailing  vortex 
system  seems  to  be  unwarranted  in  view  of  the  results  shown  in  this 
section.  Therefore  the  forces  and  moments  due  to  the  gust  penetra¬ 
tion  and  the  vehicle's  motion  are  calculated  by  using  a  modified  strip 
theory.  The  strip  theory  concept  is  a  means  cf  applying  two-dimen¬ 
sional  aerodynamic  theory  to  wings  of  finite  span.  The  more  rigorous 
lifting  surface  theories  have  largely  superseded  the  strip  theory 
approach.  However,  approximate  methods  such  as  the  strip  theorv  otill 
remain  useful  for  practical  applications  where  their  simplicity, 
flexibility,  and  economy  are  advantageous. 

The  strip  theory  method  assumes  that  the  aerodynamic  forces  or. 
each  strip  are  taken,  to  be  those  associated  with  a  two-dimensional 
wing  undergoing  the  sane  motion  as  the  particular  section  of  a  finite 
wing.  This  assumption  implies  that  tie  aerodynamic  interaction 
between  strips  is  ignored.  This  simplification  is  not  a  very  good 
approximation  of  the  three-dimensional  loading  near  the  wing  tip.- 
However,  it  is  acceptable  if  the  aspect  ratio  is  high. 

The  loads  developed  on  each  strip  can  be  calculated  from  the 
followirg  formulas  described  in  References  45-48. 

Ft 

i  / 

'  ®a  ~  4  u  c  l  a  I  __JL  v(t-r)d- 
<J  Q  d  7 


a  [w, 
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(o)i  (t) 


dr] 
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1,  -  11  P;-1-*  1  l!  9 

X  4 

*2  -  —y-2"  -  [-W  +  (c/ 2  -  Xcg)0] 

where 

mg  =  aspect  ratio  correction 

l  *  width  of  the  panel 

a  *  slope  of  the  lift  curve 

W „  *  gust  velocity 

O 

W3c  *  vehicle  velocity 

T 

xcg  ■  distance  to  c.g.  location 

$(t)  ■  Wagner  function 

•;(c)  *  Kussner  function 

The  aspect  ratio  correction  is  obtained  fior.  the  ♦  omul  a. 

m  „  _ fo_ _ 

a  1  +  (a0/::ARe)  (1+r) 

ARe  *  AR  for  symmetrical  lift  distribution 
ARe  *  AR/2  for  antisymmetrical  lift  distribution 
i  =  Glauert's  correction  factor 

The  above  integral  equations  can  be  put  imo  a  form  acre  suitable 

for  programming  by  replacing  the  integrals  by  summations . 

i  n=t  j , 

ig  =  “a  y  c  U  c  1  a  l  AWCtj) 


59 


l 

m 


mapUcJa  [W3c'o)  *(t) 

~ 

n® :  / 

I  )  •  (t-H) 

c  4 


By  summing  the  contributions  of  each  panel  the  total  lift  can  be  ob¬ 
tained.  In  a  similar  manner  the  other  forces  a--1  moments  can  be  cal¬ 
culated  . 

The  general  form  of  the  aquations  used  to  calculate  the  forces 
and  moments  acting  on  the  aircraft  were  divided  into  four  parts.  For 
example  the  Z  force  is 


Z(t)  =*  ZQ  +  Zg (t)  -I-  Z^t)  +  Z- ( t) 

where  the  sub  c,  g,  m  and  c  refer  to  the  contribution  due  to  the 
trimmed  condition,  gust  penetration,  vehicle  motion  and  control  input 
respectively.  The  components  are  defined  as 

rt  r  b/2 


Zg(t) 


hZgw(t»y)j^-  (tlty)  dy  at1 
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dy  dt. 
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Zc(t)  =  Cz j  5e(t)  QS 
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where 


2gw 


a  Q  c(v) 


kgUp 


a.  Q  c(y) 

U 


Mci> 


h  -  a-  2  c&) 
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It  was  determined  by  numerical  experimentation  that  the  forces  and 
moments  acting  on  the  aircraft  due  to  the  vehicle's  motion  could  be 
accurately  modeled  by  using  the  classical  dynamic  stability  deriva¬ 
tives.  Thus*  the  Integrals  representing  2^(6) *  Ym(t) ,  1^(1) »  etc. 
were  replaced  by  the  damping  derivatives. 

As  the  aircraft  approaches  the  vortex  core  the  induced  veloci¬ 
ties  at  some  of  the  control  points  could  become  quite  large.  The 
equations  developed  in  this  section  would  then  over  predict  the 
forces  acting  on  the  aircraft.  In  order  to  prevent  this  from 
occurring  a  ch^ck  was  incorporated  into  tl-a  aerodynamics  subroutine. 
The  check  simply  limited  the  maximum  sectional  lift  coefficient  to 
the  sectional  C^ay .  This  assumption  is  not  very  restrictive  when 
on'1  considers  the  velocities  of  the  probe  aircraft.  For  this  study 
the  probe  aircraft  had  approach  velocities  in  the  neighborhood  of 
200  I  sec.  Therefore  very  large  induced  velocities  could  be 
tolerated  before  the  angle  of  attack  approached  the  stall  region. 


Equations  of  Motion 

The  rigid  body  equations  of  motion  of  an  aircraft  having  six 
degrees  o  freedom  and  referred  to  a  rotating  axes  system  can  De 
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obtained  from  the  following  vector  equations  as  presented  by  Etkin 


bl 


F  «  m  —7—  +  mu  x  Vr 
at 

M  =  -j—  +  u  x  H 
dt 

The  rotating  axes  system  is  fixed  to  the  aircraft’s  center  of  gravity. 
The  main  advantage  of  using  a  body  axes  system  is  that  the  moments  and 
products  of  inertia  remain  constant  with  respecc  <.0  the  body  axes 
system.  The  coordinate  systems  are  shown  in  Figure  25.  Now  the 
scalar  equations  car  be  written  in  the  dimensional  nonlineari^ed  form 
as  follows: 


X(t)  -  m  g  sin  9  -  m  (U  +  QW  -  RV) 

Y(t)  *  m  g  cos  9  sin  $  *  m  (V  +  RU  -  PW) 

Z(t)  +  m  g  cos  9  cos  $  »  m  (V  +  PV  -  QU) 

L(t)  Ij^P  +  QR  (I2Z  -  Lyy) 

M(t)  »  IyyQ  +  RP  (I,,  -  IZZ) 

N(t)  =  IZZR  +  PQ  ( lyy  ’  lxx> 


The  Euler  angles  (v,  9,  $ )  can  be  obtained  bv  integrating  the  Euler 
rates  which  are  expressed  in  terms  of  the  angular  velocities  [P,  Q,  R] 
by  the  transformation 


1— •  “ 

—  *“ 

0  cos  $  -sin  i 

P 

1  sin  4>  tan  9  cos  <j>  tan  6 

Q 

0  sin  $  sec  e  cos  $  sec  9 

R 

L_ 

_  — 

The  coordinates  of  the  aircraft  relative  to  the  inertial  frame  are 
^pressed  as  a  function  of  the  Euler  angles  and  the  velocity  com¬ 
ponents  [U,  V,  Wj. 
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where  C  =  cos  S  =  sin 

The  forces  and  moments  acting  on  the  aircraft  were  divided  into 
four  components 
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Usually  the  initial  conditions  were  such  that 
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Control  -  Human  Pilot 

One  of  the  earliest  attempts  to  model  mathematically  the 
behavior  of  a  human  being  was  conducted  by  Diamantides'^  at  the 
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Goodyear  Aircraft  Corporation.  Diamantides  and  his  associates 
developed  a  pilot  model  with  the  aid  of  an  analog  computer.  This  was 
accomplished  by  postulating  various  mathematical  models  to  describe 
the  human  dynamics  and  then  by  comparing  the  results  of  the  analog 
with  human  pilot  actions.  The  test  consisted  of  a  movable  simulator 
equipped  with  a  control  column  and  display,  ihe  longitudinal  eq  ia- 
tions  of  motion  for  a  fighter  aircraft  were  used  to  test  the  pilot's 
reaction  to  pitch  disturbances.  The  human  pilot  and  the  analog  were 
run  in  parallel.  This  allowed  easy  evaluation  of  the  postulated  human 
dynamics.  The  test  engineers  adjusted  the  various  parameters  in  the 
analog  model  until  the  pilot  was  unable  to  detect  that  the  analog  was 
flying  the  simulator.  Figure  30  shows  a  block  diagram  representation 
of  Diamantides'  pilot  model.  This  figure  indicates  that  upon  per¬ 
ception  of  the  stimulus  the  human  controller  performs  two  linear 
operations.  First  he  makes  a  mental  computation  of  the  stimulus. 

That  is,  position,  rate  and  acceleration  are  sensed  and  weighed, 
then  a  decision  to  act  occurs.  Secondly,  there  art  certain  physical 
limitations  on  the  corrective  hand  motions  which  are  caused  by  the 
neuromuscular  or  motor  feedback  loop.  Another  characteristic  of 
the  human  operator  is  that  he  normally  does  not  react  to  a  stimulus 
until  it  is  above  a  certain  threshold.  Once  the  magnitude  of  the 
stimulus  is  higher  than  the  threshold  value  the  information  is  sensed 
and  a  mental  computation  is  performed.  The  brain  then  commands  the 
appropriate  muscle  reaction.  This  all  takes  place  in  a  fraction  of 
a  second.  The  time  lag  is  called  the  "reaction  time".  It  has  been 
found  that  the  reaction  time  is  almost  constant  in  normal  people. 
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There  are  two  types  of  control  possible.  They  are  pursuit  and 
compensatory.  In  pursuit  tracking  past  experience  provides  the 
tracker  with  information  about  what  to  expect  in  a  future  input.  In 
compensatory  tracking  this  is  not  possible^.  Figure  31  is  a  block 
diagram  representation  of  both  pursuit  and  compensatory  tracking. 

In  studying  the  dynamic  behavior  of  an  aircraft  penetrating  a  v  tex 
system  it  is  quite  apparent  that  the  compensatory  control  is  tht- 
oniy  type  which  is  applicable. 

The  «ilot  model  for  compensatory  control  case  can  be  expressed 
in  terms  of  the  following  transfer  function. 


V„<s) 


±2.  e 


-ts 


(1  +  tls) 


(1  +  Tns)(1  +  TjS) 


kp  =  pilot  gain 
e  Tb  »  reaction  time  delay 
-  lead  time  constant 
“  neuromuscular  lag  time  constant 
Tj  *  lag  time  constant 


The  above  expression  can  be  put  into  a  mere  convenient  form  by 
assuming  that  the  neuromuscular  lag  can  be  accounted  for  by  modifying 


the  reaction  time  delay. 

(Ii£_±l>  -~TS 

CpS  +  1) 


yp(s)  “  *p  (77 ”•  e 


The  parameters  included  in  these  expressions  should  not  be  inter¬ 
preted  as  constants.  In  addition  to  the  transfer  function,  adjust¬ 
ment  rules  are  required  to  specify  the  pilot  model  parameters. 

These  rules  have  been  determined  by  numerous  investigators . ^ 


It.  has  been  found  that  the  gain  kp  is  easily  and  rapidly  varied 
by  the  pilot  to  an  appropriate  value.  The  equalization  tern,  that  is 
(t^s+1)  /  (ipS+l),  is  varied  within  certain  limits  so  tha'-  the  system 
has  a  good  frequency  response. 

The  adjustment  rules  can  be  stated  as  follows: 

1.  The  human  adapts  so  that  his  gain  ana  equalization  charac¬ 
teristics  are  appropriate  for  a  stable  system. 

2.  The  human  adapts  so  that  the  form  of  his  equalization  char¬ 
acteristics  is  appropriate  to  good  low-frequency  closed-loop  response 

to  the  forcing  furction. 

The  adjustment  rules  were  determined  by  analyzing  experimental  data 
obtained  from  tests  on  human  operators.  To  put  the  cdjustmen,.  rules 
in  a  more  practical  format  we  can  use  the  following: 

1.  The  open-loop  phase  margin.  The  pilot  adapts  3n  equaliza¬ 
tion  to  obtain  50  to  110  degrees  of  phase  margin.  It  has  been  found 
that  an  open-loop  phase  margin  in  this  range  yields  a  closed  loop 
system  with  a  satisfactory  time  uomain  performance, 

2.  Cpen-loop  crossover  frequency.  For  good  high  frequency 
characteristics  the  crossover  frequency  should  be  at  least  1  rad/sec. 

I'sing  these  rules  the  constants  tl  and  tj  can  be  determined  fioc 

eithe’-  a  Bode  or  Polar  plot  of  the  open  loop  transfer  function 

fY..  Y  1,  wnere  Y„  is  the  transfer  function  of  the  aircraft.  Another 

pa  » 

technique  vhicn  can  be  employed  to  determine  the  pilot  parameters 
is  the  root  locus  method. 

Both  single  ana  multi  axis  control  were  incorporated  into  the 
program.  Based  on  the  work  in  References  A?,  AS,  and  53-59  it  was 


CT\ 


69 


decided  that  the  pitch  and  roil  axis  controllers  could  be  modeled 
accurately  by  the  simple  transfer  function 


Pa  K  kpJTL4s+1>e"TS 

r  *  v  v 

(Roll) 

PO  *  kne(''L.-s+1)e~  S 

(Pitch) 

Figure  32  sr.ows  a  block  diagram  representation  of  tne  pilot-vehicle 
system.  Tne  transfer  function  can  be  written  in  the  time  domain  as 


MtvO  -  kp9  ee(t)  +  kPe  tL(, 


d'J 


dt 


anc 


dis 


+  kpt  'ti  jr 


where 


*■  <fc  -  $■  «  -  1  for  wings  level  control 

Se  *  ?c  -  0. 

Appendix  A  gives  the  details  on  the  selection  of  the  pilot  parameters. 
Figure  33  illustrates  the  block  diagram  representation  of  the  roll 
and  heading  control  task.  Again  the  pilot  model  was  assumed  to  have 
the  same  form  as  those  previously  described. 

The  pilot  oodels  used  in  this  study  should  give  a  qualitative 
assessment  of  the  pilot's  ability  to  control  the  upset  due  to  a 
vortex  encounter. 


Numerical  Analysis 

The  twelve  first  order  differential  equations  needed  to  describe 
the  aircraft's  motion  were  solved  numerically  by  using  a  fourth  or-’er 
Runge-Kucta  technique.  The  Runge-Kutta  method  for  a  system  of 
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differential  equations  can  be  found  in  any  text  on  numerical  analysis. 
For  completeness  the  general  form  of  the  Runge-F.utta  algorithm  is 
presented  below  in  its  functional  form.^ 


Pk  '  h  fk(Xi,  yii,  y2>1»  •  *  •  •  yri>i> 

qk  =  h  fk(Xi  +  h/2,  yx>i  +  Y'  *  '  *  *  *  yn,i  +  F} 

rk  -  h  fk(Xi  +  h/2>  ylfl  +  *  *  *  '  .  yn,i  +  F} 

sk  -  h  fk(Xi  +  h,  ylfl  +  xL  ,  •  •  •  •  ,  yQ>i  +  rn) 

yk,  i+i  •  ykfi  +  X/6  (Pk  +  2%  +  2rtc 4  sk> 


where  k  «=  1,  2,  3  •  •  •  .  h 

The  equations  of  motion  rearranged  to  a  form  suitable  for  integration 
are  as  follows: 


U  »  (X(t)  -  m  g  sin  6)/m  -  QW  +  RV 

V  «  (Y(t)  +  m  g  cos  0  sin  o)/m  -  RU  +  PW 

W  «  (Z(t)  +  m  g  cos  0  cos  $)/m  -  PV  +•  QU 


P  -  (L(t)  -  QR  (Izz  -  lyy))/!^ 

Q  -  (M(t)  -  HP  (1„  -  Izz>)/V, 

R  «  (N(t)  -  PQ  (Iyy  -  lx*))/!** 

0  -  Q  cos  $  -  R  sin  $ 

$  “  •  +  0  sin  ?  tan  6  +  R  cos  $  tar.  0 


^  *  Q  sin  t  sec  0  +  R  cos  $  sac.  P 

X  -  U  C  C  +  V(S  SnC  -  ChS  )  +  W(C  S^C  +55,1 
0  V  v  8  V  0  *■  >0ii)  51; 

Y  -  U  C.S,  +  V(S.S.S,  +  C„Ct)  +  W(C.S„S.  -  sc,) 

2  -  U  s0  +  V(s^ce)  +  w<ey;.) 

where  C  “  cos  S  *  sin 
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SECTION  IV 


RESULTS 


Simplified  Analysis 

Before  discussing  the  results  of  the  computer  simulation  it 
would  be  Instructive  to  examine  the  results  from  a  simplified  analysis. 
Consider  the  situation  of  an  aircraft  suddenly  penetrating  along  the 
axis  of  the  vortex  core.  A  measure  of  the  vortex  hazard  would  be  a 
comparison  of  the  vortex  induced  rolling  moment  to  that  available  due 
to  maximum  aileron  deflection.  This  type  cf  analysis  has  been  dis¬ 
cussed  previously  in  Chapter  1.  Figure  34  taken  from  reference  12 
indicated  the  severity  cf  such  a  vortex  encounter.  The  vortex  gene¬ 
rated  by  the  leading  wing  is  approximately  that  of  a  DC-9  during 
approach.  The  separation  distance  and  separation  time  were  9000  feet 
and  one  minute  respectively.  Examination  of  this  figure  reveals  that 
penetrating  aircraft  would  experience  induced  rolling  moments  exceed¬ 
ing  their  roll  capability.  Figure  35  also  taken  from  reference  12 
shows  the  effect  of  the  lateral  displacement  on  the  vortex  induced 
rolling  moment.  Note  that  the  induced  rolling  moment  reverses  sign 
as  the  vo i ten  ccves  toward  the  wing  tip. 

Using  simple  strip  theory  the  induced  rolling  moment  was 
calculated  for  aircraft  of  varying  sizes.  Table  1  contains  the 
characteristics  of  the  penetrating  aircraft  while  Table  2  contains 
the  information  pertaining  to  the  generating  aircraft.  The  results 
of  theca  calculations  are  presented  in  Figures  36  -  38.  Figure  36 
shows  that  a  light  single  engine  aircraft  would  experience  momentary 
loss  of  control  for  separation  distances  of  up  to  9  miles  behind  a 
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DC-6.  Also  the  light  twin  would  be  cut  of  control  for  penetrations 
under  2.3  miles  separation.  Figure  37  shows  that  both  the  lighc 
single  and  twin  engine  aircraft  are  momentarily  out  of  control  behind 
a  DC-9.  Note  that  for  the  light  twin  to  be  under  control  the  separa¬ 
tion  distance  must  be  in  excess  of  3  miles.  Now  in  Figure  38  we  see 
that  even  an  aircraft  the  size  of  a  DC-9  would  be  consid 'red  to  be  ou 
of  control  for  separation  distances  of  up  to  13.5  miles  behind  a  747. 
Figure  39  taken  from  ?  repor-  by  Robinson  and  Larson*^  illustrates 
that  the  strip  theory  apprc^h  yields  an  effective  upper  bound  to  the 
induced  rolling  taoment.  The  scatter  in  the  flight  test  data  is  due 
largely  to  the  difficulty  of  penetrating  into  the  vortex  core.  This 
will  be  discussed  in  more  detail  later  in  this  section.  An  interest¬ 
ing  result  of  the  flight  test  program  conducted  by  NASA  was  that  the 
pilots  felt  that  they  had  control  of  the  aircraft  at  separation 
distances  vhe.re  the  ratio  of  C^/C^*  5^^  was  greater  than  one.  This 
is  due  most  likely  to  two  factors.  First  as  the  aircraft  begins  to 
roll  the  natural  roil  damping  due  to  the  wing  aids  the  pilot  in 
regaining  control  of  the  aircraft.  Secondly  and  probably  the  most 
important  is  that  the  aircraft  is  rapidly  expelled  from  the  vortex 
system  thus  giving  the  pilots  the  false  Impression  that  they  were  in 
c  -dDlete  control  of  the  aircraft. 

Along-Track  Penetration  With  Roll  Control 

Along-track  penetration  is  most  likely  to  occur  during  final 
approach  or  shortly  after  take-off.  In  either  case  the  following 
aircraft  can  ill-afford  to  ha^e  any  large  excursions  from  it--;  normal 


Figure  37.  Induced  Rolling  Moment  Due  to  the  Wake  of  a  DC- 
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flight  path.  For  the  separation  titles  considered  in  this  study  the 
normal  downward  motion  of  the  vortex  system  \  ould  place  the  vortices 
well  below  the  flight  path  of  the  following  ircraft.  However,  several 
investigators  have  shown  that  under  certain  tmospheric  conditions  the 
downward  motion  of  the  vortex  system  can  be  re  arded.  In  fact  one 
investigator  has  suggested  that  as  the  vortex  system  approaches  an 
inversion  layer  the  vortices  behave  as  though  they  were  approaching 
the  ground.  Thus  the  vortices  would  cease  their  downward  movement 
and  begin  to  spread  laterally.  Another  possibility  of  encountering  a 
vortex  wake  in  the  terminal  area  would  he  during  the  execution  of  a 
missed  approach  by  the  lead  aircraft.  The  wake  of  the  leaa  aircraft 
would  be  laid  down  at  an  altitude  well  above  that,  of  the  following 
aircraft  ard  thus  could  descend  into  the  flight  path  of  the  trailing 
aircraft. 

Now  consider  some  typical  results  from  the  six  degree  of 
freedom  computer  simu>atlon.  First  let  us  examine  the  response  of  an 
air '.raft  where  the  pilot  controls  the  roll  attitude.  That  is,  the 
mathematical  pilot  model  attempts  to  maintain  a  wings  level  orienta¬ 
tion.  The  char cu.t eristics  of  the  penetrating  and  generating  aircraft 
are  listed  In  Table  III.  Figure  40  illustrates  the  penetrating  angles 
under  consideration.,  is  the  angle  between  the  penetrating  air¬ 
craft’s  flight  path  and  the  plane  of  the  vortices.  Whereas  ''p  is 
defined  ss  the  angle  between  the  flight  path  and  the  vorten  axis  when 
viewed  from  above.  For  the  cases  considered  here  the  separation 
distance  was  27,000  feet  or  approximately  2  minutes  behind  the 
generating  aircraft.  Alao  it  should  be  pointed  out  that  the  flight 
path  of  the  penetrating  aircraft  was  selected  so  that  it  would 
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intersect  the  vortex  lore  if  the  aircraft  were  not  influenced  by  the 
vortex  system. 

For  an  aircraft  descending  into  the  vo-  -ex  system  at  the  =»;»gle 
Rp  v  2.5J,  Figures  41  -  43  illustrate  the  «ffe *t  of  "He  oblique  pene¬ 
tration  angle  Vp  on  the  response  of  the  penetrating  aircraft.  Notice 
that  the  aircraft  descending  into  the  vortex  with  pp  =  3*  .oils  away 
from  the  vortex  system.  The  aircraft  is  seen  to  roll  slightly  and 
climbs  away  from  the  vortex.  The  aircraft  climbs  since  it  is  flying 
into  the  upwash  region  of  the  port  vortex.  Upon  increasing  the  pene¬ 
tration  angSe  to  ^p  -  w.»  see  a  significant  difference  in  the  air¬ 

craft's  response.  Figure  42  shows  the  aircraft  rolling  sharply  to  the 
right  and  at  the  same  time  its  rate  of  descent  is  increased  as  the. 
aircraft  passes  into  the  strong  downwash  region  between  the  /ort_~es. 
Figure  43  illustrates  the  effect  of  including  the  pilot  in  trw.  simula¬ 
tion.  Without  pilot  control  the  aircraft  rolls  to  an  angle  [.-eater 
than  100  degrees  whereas  with  cuntrcl  input  the  aircraft  rolls 
approximately  50°  before  roll  recovery  takes  place.  The  time  history 
plot  of  the  yaw  angle  reveals  the  influence  that  the  vertical  tail 
plays  tr  the  aircraft's  dynamic  behavior,.  As  the  aircraft  enters  the 
vortex  field  the  aircraft  is  seen  to  roll  in  a  counterclockwise 
direction.  This  can  be  explained  by  referring  back  to  Figure  35 
which  illustrates  the  affect  of  lateral  displacement  on  the  induced 
rolling  moment  coefficient.  Note  that  the  iaco'ced  rolling  moment 
coefficient  reverses  sign  as  the  vortex  moves  toward  the  wing  tip. 

Thus  as  the  aircraft  approaches  the  port  vortex  it  initially  rolls  in 
the  opposite  direction  with  respect  to  the  vortex  field.  The  pilot's 
reaction  to  the  disturbance  would  be  to  apply  aileron  control  to  roll 
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the  aircrafc  in  a  clockwise  direction.  Now  as  the  aircraft  naar°  the 
vortex  the  aircraft  Is  rolled  sharply  in  the  clockwise  direction.  Thus 
we  see  that  the  pilot's  initial  reaction  would  momentarily  be  out  of 
phase  with  the  roll  disturbance.  Also  during  the  vortex  encounter  the 
vertical  tail  experiences  an  induced  velocity  irom  the  left  which 
causes  the  aircraft  to  yaw  to  the  left  by  approximately  12°.  Vs  he 
aircraft  passes  into  the  field  of  influence  of  the  starboard  vortt 
the  upsetting  roll  moment  tends  to  aid  tr.e  pilot  in  regaining  a  wings 
level  attitude.  The  yaw  excursion  to  the  right  is  due  primarily  to 
the  roll  orientation. 

Figures  44  -  49  show  th*>  effect  of  varying  the  angle  6p.  For 
these  calculations  ijip  is  held  constant.  Here  the  aircraft  is 
climbing  into  the  vortex  system.  The  angle  0p  was  varied  from  3° 
to  9°.  These  figures  again  indicate  the  influence  of  the  vertical 
tail  on  the  aircraft’s  response.  As  the  aircraft  climbs  iotc  the 
velocity  field  of  the  port  vortex  the  vertical  tail  experiences  a  side 
velocity  from  the  right.  This  gives  rise  to  a  positive  yawing  moment 
which  yaws  the  airplane  to  the  right.  As  the  aircraft  passes  through 
the  plane  of  the  vortices  its  rate  of  climb  diminishes  and  the  air¬ 
craft  roils  in  a  clockwise  direction.  The  aircraft  in  Figure  44  is 
approximately  60  feet  below  and  80  feet  to  the  left  of  the  port 
vortex.  However,  after  8  seconds  the  aircraft  is  located  49  feet 
benea'h  the  vortex  plane.  If  the  airplane  had  not  encountered  the 
vortex  system  the  airplane  would  have  been  approximately  40  feet 
above  the  plane  of  the  vortices.  Similar  resuits  are  shown  in 
Figures  44  -  49.  Based  on  the  results  presented  here  several  conclu¬ 
sions  can  be  made  concerning  the  dynamic  behavior  of  the  following 
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aircraft.  The  calculations  shoved  that  even  with  pilot  control 
large  excursions  in  roll,  yaw,  and  altitude  were  experienced  b>  the 
penetrating  aircraft.  Also,  the  results  revealed  the  significant 
influence  that  the  vertical  tail  makes  to  the  aircraft's  dynamic 
response.  Finally  for  the  cases  considered  in  this  section  the 
following  aircraft  was  in  close  proximity  to  the  ground  and  thus  the 
excursions  shown  here  are  completely  unacceptable  from  a  flight  safety 
standpoint. 

Additional  results  of  the  along  track  penetration  study  are  show, 
in  Figures  50  -  53.  The  probe  aircraft  selected  were  a  Jet  Star  and 
DC -9.  These  aircraft  were  flown  into  the  wake  of  either  a  DC-9,  727 
or  a  DC-10.  Figure  50  shews  the  maximum  roll  excursion  and  roll  rate 
of  the  Jet  Star  as  a  function  of  separation  distance.  The  figure 
indicates  that  there  is  very  little  difference  in  the  aircraft's 
response  at  separation  distances  from  2  to  8  miles.  This  can  he 
explained  by  examining  Figure  51  which  is  a  plot  of  the  maximum 
induced  rolling  moment  coefficient  as  a  function  of  separation  dis¬ 
tance.  The  equation  used  in  computing  this  curve  was  taken  from 
reference  12.  Note  that  the  rolling  moment  coefficient  decreases 
slowly  with  increasing  separation  distance.  This  is  due  to  the 
method  used  to  predict  the  vortex  velocity  distribution.  Recall  that 
McCormick's  equations  predict  small  core  radii  and  large  maximum 
tangential  velocities.  The  lower  portion  of  Figure  51  3hows  the 
velocity  distributions  used  in  computing  the  roll  coefficient.  At 
distances  of  several  core  radii  along  the  span  the  induced  velocities 
tend  to  coalesce.  This  is  the  region  which  contributes  the  major 
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portion  of  the  induced  rolling  moment.  Therefore  it  follows  that  the 
roil  coefficient  will  decrease  slowly  with  separation  distance. 

As  mer.  tic  ted  in  the  introduction  various  techniques  have  been 
proposed  to  modify  :he  vortex  structure.  Some  of  the  concepts  which 
have  been  considered  include  spoiler  deflection  and  mass  injection 
into  the  vortex  core.  Both  of  these  techniques  reduce  the  maximum 
tangential  velocity  and  increase  the  core  radius.  However,  at  a  dis¬ 
tance  of  several  core  radii  from  »he  vortex  center  the  velocity 
distributions  are  the  same  as  the  unaltered  vortex.  This  can  be  seen 
in  Figure  11.  This  figure  presents  Hardman's  data  and  illustrates 
the  effect  that  blowing  has  on  the  vortex  velocit"  distribution.  Based 
on  the  results  shown  in  Figure  51  it  can  be  concluded  that  the  vortex 
abatement  devices  will  not  eliminate  the  vortex  hazard  unless  the 
device  hastens  the  vertex  break-up. 

In  the  literature  it  has  been  suggested  that  the  vortex  hazard 

can  be  assessed  by  the  ratio  of  the  maximum  induced  acceleration  to 

the  maximum  acceleration  obtained  with  full  aileron  deflection.  If 

the  ratio  1 P ] / Pa  is  greater  than  1  the  aircraft  is  assumed  to  be 

'  "max  ° 

out  of  control.  This  ratio  is  of  course  the  same  as  the  ratio 
Ci/Ci.  °a  used  previously.  Figure  52  shows  that  the  Jet  Star  flying 
into  the  wake  of  either  a  DC-9  or  / 27  would  be  uncontrollable  for  the 
separation  distances  shown.  The  roll  excursion  of  the  Jet  Star  in  the 
wake  of  the  DC -10  or  727  exceeded  200  degrees  ror  the  same  separation 
distances.  Now  Figure  53  shows  the  same  type  of  results  for  a  DC-9 
in  the  wake  of  either  a  DC-9,  727  DC-10.  The  maximum  roll  angles 

exceed  40°  for  penetration  into  the  wake  of  the  727  and  DC-10. 
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Whereas  the  maximum  roll  angle  in  the  wake  of  the  DC-9  is  only  20°. 

The  ratio  Iplcal/Pfi  shown  in  Figure  56  indicates  that  the  DC -9 
max 

is  nomentarily  out  of  control  when  penetrating  the  wake  of  a  727,  DC -9 
and  DC -10. 

In  an  effort  to  assess  the  validity  of  the  computer  simulation 
the  results  were  compared  with  flight  test  data  obtained  from  Refer¬ 
ence  62.  These  data  were  obtained  by  the  NASA  Flight  Research  Center. 
The  probe  aircraft  were  positioned  by  radar  to  a  specified  distance 
behind  the  generating  aircraft.  The  probe  aircraft  were  then  flown 
into  the  vortex  wake  for  2  to  3  minutes  in  order  to  record  sufficient 
response  data.  The  vortex  system  of  the  generating  aircraft  were 
marked  by  entrainment  of  the  smoke  generated  from  the  aircraft's 
engines.  However,  it  was  noted  that  there  was  a  lack  of  smoke  to 
clearly  mark  the  wake  at  large  separation  distances. 

The  computed  results  for  a  DC-9  penetrating  the  wake  of  a 
Contair  990  are  presented  ir.  Figure  55  along  with  the  flight  test  data. 
Two  calculated  curves  are  shown  in  this  figure.  The  solid  curve  is 
based  on  calculations  including  only  roll  control  and  the  dashed  curve 
for  both  roll  and  heading  control.  The  solid  curve  is  low  in  compari¬ 
son  to  the  flight  test  data  for  separation  distances  under  4  miles. 

In  the  5  to  S  mile  range  the  computed  results  are  high.  However,  the 
flight  test  data  exhibits  considerable  scatter.  This  is  due  largely 
to  the  difficulty  of  penetrating  the  vt rtex  system.  Recall  Figure  35 
which  illustrates  the  effect  of  lateral  ^paration  on  the  vortex 
induced  rolling  moment.  Inis  figure  indicates  that  if  the  aircraft 
missed  the  vortex  core  by  a  relatively  small  distance  the  roll  ex¬ 
cursion  would  he  quite  different.  Since  the  pilots  had  difficulty 
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la  seeing  the  vortex  wake  at  large  separation  distances  one  would 
expect  to  have  scatter  in  the  flight  test  data.  Another  factor  con¬ 
tributing  to  the  scatter  is  the  break  up  of  the  vortex  system. 

The  influence  of  trying  to  maintain  heading  is  illustrated  by 
the  dashed  curve.  By  maintaining  heading  the  pilot  increases  the  roll 
upset.  As  previously  shown  the  aircraft  is  rapidly  expelled  from  the 
vortex  system  and  thus  by  trying  to  hold  his  heading  the  pilot  keeps 
the  aircraft  in  the  influence  of  the  vortices  for  a  longer  time.  This 
then  gives  rise  to  larger  roll  excursions.  Based  on  this  comparison 
it  is  felt  that  the  results  of  the  computer  simulation  are  representa¬ 
tive  of  the  actual  aircraft-vortex  encounters. 

Figure  56  is  a  time  history  plot  of  the  roll  and  yaw  angle.  The 
pilot  model  is  successful  in  maintaining  the  aircraft's  heading.  In 
this  case  the  roll  excursion  is  essentially  the  same.  The  difference 
between  this  and  the  preceding  case  is  the  penetration  angle  -^p.  In 
the  second  example  the  pilot  is  attempting  to  maintain  an  oblique 
angle  to  the  vertex  axis.  However,  in  the  results  shown  in  Figure  55 
the  pilot  is  attempting  to  fly  right  along  the  vortex  axis. 

Numerous  cases  of  cross  track  penetration  were  calculated.  The 
results  shown  in  Figures  57  and  58  are  representative  of  the  cases 
studied.  The  normal  load  factors  without  control  iixput  compare 

A  5  6 

favorably  with  the  results  obtained  in  previous  investigations.  *  * 

However,  when  pilot  control  was  added  the  normal  load  factor  was 
increased  slightly  which  is  quite  different  from  the  results  shewn  in 
igure  3.  The  results  of  this  study  indicate  that  structural  failure 
is  not  a  serious  threat  to  tha  penetrating  aircraft.  Figure  58  shows 
the  angle  of  attack  time  history.  Notice  that  the  wing  is  stalled 


Business  aircraft  in  the  wake 
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for  a  fraction 
to  develop  any 


of  a  second;  however,  the  aircraft  does  not  have  time 
appreciable  vertical  velocity. 
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SECTION  V 
CONCLUSIONS 

Based  on  the  results  presented  in  the  previous  section  the 
following  conclusions  have  been  reached  concerning  the  dynamic  behavior 
of  an  aircraft  penetrating  a  vortex  wake. 

1.  Encountering  a  vortex  can  produce  rolling  moments  on  the 
aircraft  in  either  direction. 

2.  The  pilot's  reaction  to  the  roll  disturbance  can  be  momen¬ 
tarily  out  of  phase  with  the  upsetting  roll  moment. 

3.  An  aircraft  descending  into  a  vortex  can  experience  com¬ 
pletely  different  types  of  response  depending  upon  the 
penetration  angle  yp.  For  very  shallow  angles  the  aircraft 
is  rolled  away  from  the  vortex  system.  As  increases  the 
aircraft  is  rapidly  relied  to  a  very  large  roll  orientation. 

4.  Aircraft  climbing  into  the  vortex  system  will  experience 
large  roll,  yaw  and  altitude  excursions. 

5.  Reduction  of  the  maximum  vortex  tangential  velocity  does  not 

significantly  reduce  the  hazard  of  a  .-  —ex  encounter. 

6.  The  calculated  results  agree  favorably  with  flight  test 
data  for  up  to  5  miles.  At  distances  greater  than  5  miles 
the  computed  results  were  too  high. 

7.  The  pilot  can  Increase  tne  roll  upset  by  trying  to  maintain 
his  heading.  In  so  doing  the  aircraft  remains  in  the 
vortex  field  for  a  longer  time  resulting  in  a  larger  roll 
excursion. 

8.  In  the  transverse  penetration  the  pilot's  control  input  did 
not  appreciably  increase  the  load  factor  as  shovi  iu  previous 
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studies.  For  the  aircraft  used  in  this  study  the  results 
indicate  that  structural  failure  is  not  a  serious  threat 
to  the  penetrating  aircraft 

9.  The  results  show  that  corporate  and  light  jet  transports 
cannot  mix  safely  with  the  heavier  jet  tran'-ports  at  dis¬ 
tances  less  than  at  least  8  miles. 

10.  Admittedly  the  probability  of  encountering  a  vortex  wake  “s 
still  quite  smallj  however,  if  the  separation  distances  are 
reduced  the  possibility  of  an  encounter  will  definitely 
increase.  The  results  of  this  study  indicate  that  reduced 
separation  distances  are  unacceptable  from  a  flight  safety 
viewpoint. 
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Appendix  A 

Selection  of  Pilot  Parameters 

The  pilot  parameters  we’-e  selected  on  the  basis  of  the  rules  pre¬ 
sented  in  Chapter  III.  In  order  to  apply  these  rules  the  open  Loop 
transfer  function  for  the  desired  control  task  must  oe  develfpe 
The  equations  used  in  selecting  the  pilot  parameters  are  present-- ' 
for  the  following  control  tasks: 

1.  Pure  Poll 

2.  Roll  and  Heading  Control 

3.  Pure  Pitch 


Roll  Control 

The  equation  governing  the  rolling  motion  of  an  airplane  can  b 
expressed  in  terms  of  the  Laplace  transform  variable  as 
(s2-LpSH  -  L6a6a 

or  rewritten  to  form  the  roll  angle/aileron  deflection  transfer  func¬ 
tion 

*£?)...  -  _ 

6a(s)  s(s-Lp) 

Previously  it  was  stated  that  the  pilot  could  ba  modeled  by  a  gain, 
simple  lead  term,  and  a  time  delay  which  in  mathematical  notation 


f unctions . 
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It  has  been  found  by  experimentation  that  pilots  adopt  an  equaliza¬ 
tion  term  so  that  they  cancel  the  roll  mode.  That  is  the  pilots  lead 
term  is  fixed  by  the  roll  root  TR  *  1/Lp.  The  pilot's  gain  Kp^  was 
calculated  from  an  expression  taken  from  Reference  b3.  For  a  well 
damped  response  Kp^  was  related  to  the  Lp  and  L$a  as  follows: 

Kn 


fl 


2U 

•"a 


Heading  and  Roll  Contr^ 

Assuming  the  roll  axis  to  be  handled  as  explained  above  the 
heading  transfer  function  can  be  approximated  by  the  following  ex¬ 
pression. 


N6(s+1/tVi) 

T  “  s(s2+2? ,uds+ud2) 

where 

Yd 

1/T  -  (Yv - No) 

No 

2»dud  “  -<Yv  +  V 

wd2  *  NS 

Now  the  open  loop  transfer  function  can  be  written. 

v_  "L.y  N’5(s+1/t.4,1)(s+1/Tl^) 

-e  s(si+2;dij(js-hjd2) 
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The  rules  stated  in  Chapter  III  call  for  a  selection  of  pilot  lead 
and  gain  so  that  the  phase  margin  lies  between  50  and  100  degrees. 
The  lead  term  can  therefore  be  adjusted  to  satisfy  the  phase  margin 
requirement  while  the  gain  can  be  varied  to  ensure  that  the  open 
loop  crossover  frequency  falls  in  the  vicinity  of  1.0  radians  per 
second.  Figure  59  shows  a  block  diagram  sketch  of  the  pilot-vc'  icle 
system  as  wall  as  a  Bode  plot  of  the  open  loop  transfer  function. 

Pitch  Control 

The  pitch  angle/ elevator  deflection  transfer  cunction  can  be 
obtained  by  solving  the  longitudinal  set  of  equations. 

(s-Xju  -  (V»  +  X^w  -  (XqS  -  g)e  -  X55 
-ZyU  +  (s-V  "  Zw)*  "  (uo  +  Zq)s0  »*  Z66 
-MyU  -  (MyS  +  My)w  +  (s2-  MqS) 6  -  M$6 

solving  for  6(s)/5(s)  yields 

e(8)  _  A9S2  +  B9s  -f  Ce 
6(a)  ”  Along 

where 

Along  *  [s^CduOpS+Up2]  UZ+lUuOgpS-h^gp2] 

The  terms  A9,  Bq,  Cq,  (Cu)p  etc.  are  related  to  the  aircraft's 
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stability  derivatives.  Aehkenas  and  McRuer  have  shown  that  the 
following  approximations  can  be  made  for  the  numerator 

„  1 

Ac.s ^  +  Bgs  +  C.%  -  As(s+^— )  (s+l/ra0) 

0  t  r.  -  2 


*  N6  (S+i/TW)<S+l/Tfcl(,) 

S(SJ4-25dwdsWr) 


Figure  59.  Bode  Plot  of  Heading  Transfer  Function 
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where 

A$  *  M<s 

zu 

1/Te1  -  -Xu  +  x*  ^ 

1/T02  - 

and  for  the  denominator 


2  6  usp  -  -  (Z„  +  Mq  +  V 

“sp 


2  {  B 
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H^Xq-g) 

Zw  Hq-Ma 


(MWZU  -  MyZy) 


^  Mq-Ma 


Again  the  open  loop  transfer  function  i9  obtained  by  multiplying  the 
aircraft  transfer  function  by  the  pilot  transfer  function. 
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[YpJ 
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ilfl  h>Q  TL9(8+l/TL6)(s+l/Tei)(s+l/T92)e  1 

0e  (82+2(6w)p8  +  uip2)(s2+  2(6u,')Sp  8  +  uSp2) 

Figure  60  shews  a  block  diagram  sketch  and  Bode  plot  of  the  open 
loop  transfer  function.  Note  that  the  crossover  frequency  is  approxi¬ 
mately  4.5  radians /second.  This  value  was  selected  based  on  the 
information  contained  in  a  paper  by  McRuer,  Graham  and  Krendel^. 

They  found  that  for  good  pitch  attitude  control  a  crossover  frequency 
of  4.5  radians /second  is  required  for  an  effective  time  delay  of 


.25  seconds. 


